Ferromagnetism in repulsive Fermi gases: 
upper branch of Feshbach resonance versus hard spheres 
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We use quantum Monte Carlo, including backflow corrections, to investigate a two-component 
Fermi gas on the upper branch of a Feshbach resonance and contrast it with the hard sphere gas. 
We find that, in both cases, the Fermi liquid becomes unstable to ferromagnetism at a ba smaller 
than the mean field result, where is the Fermi wavevector and a the scattering length. Even 
though the total energies E(kpa) are similar in the two cases, their pair correlations and kinetic 
energies are completely different, reflecting the underlying potentials. We discuss the extent to 
which our calculations shed light on recent experiments. 

PACS numbers: 67.85.-d, 37.10.Jk, 71.27.+a 



Introduction: Ultracold atomic gases are emerging 
as a unique laboratory for testing quantum many-body 
Hamiltonians. A problem of fundamental importance is 
the ground state of two species of fermions interacting 
via repulsive interactions. The attractive case is now 
well-understood and shows the BCS-BEC crossover [lj 
in the superfluid ground state. The broken symmetry is 
already apparent within BCS mean field theory (MFT) 
with an arbitrarily small attraction leading to a paired 
superfluid. In contrast, we know much less about the re- 
pulsive case. The Landau Fermi liquid, known to exist at 
weak repulsion [2j, can become unstable only beyond a 
critical value of the interaction^. Thus the phase tran- 
sition is not a weak coupling problem, and the validity 
of MFT in the repulsive case is questionable. 

An exciting new development is a recent experiment [J] 
which has been interpreted as evidence for a ferromag- 
netic instability d, U0| in a "repulsive" Fermi gas of 
e Li atoms. A crucial point is that the interactions be- 
tween the atoms are quite different from the textbook 
problem of hard-sphere interactions. In the experiment, 
the atoms are on the upper branch of a Feshbach reso- 
nance with a positive s-wave scattering length a. The 
two-body ground state then is a molecule of size a. But 
in the upper branch, where the wave function is made 
up from scattering states, the atoms feel an effective re- 
pulsion characterized by a > 0, despite the fact that the 
underlying potential is attractive. 

The main question we examine in this Letter is the ex- 
tent to which the many-body physics in the upper branch 
is similar to, or different from, that of a purely repul- 
sive Fermi gas. We use quantum Monte Carlo (QMC) to 
compute the energy, chemical potential and pair distri- 
bution function of the two systems - upper branch and 
repulsive - to understand the instability of the Fermi liq- 
uid to ferromagnetism. We believe that such a study of 
equilibrium properties is necessary, before one addresses 
non-equilibrium questions in the upper branch. 

Before describing our results, we emphasize important 
ways in which our work differs from previous studies, 



which focus on MFT of purely repulsive interactions. 
First, we carefully discuss what it means for a many- 
body wavefunction to be on the upper branch, which is 
essential to describe the experiments. Second, it is cru- 
cial to use QMC for this strong coupling problem. For 
instance, QMC calculations [8J for the electron gas show 
that ferromagnetism sets in at a critical density nearly 
3 orders of magnitude smaller than that predicted by 
Hartree-Fock MFT. Finally, we include backflow correc- 
tions, which can have a nontrivial effect on the nodes 
of the many-body wavefunction, and thus on the ground 
state energy [9j . Not including backflow may lead to spu- 
rious ferromagnetic instabilities in normal 3 He [lOj . 

Our main results are that we find Ferromagnetic (FM) 
instabilities in both the upper branch and the hard sphere 
Fermi gas. For small kpa > 0, with kp the Fermi 
wavevector and a the s-wave scattering length, both sys- 
tems are Landau Fermi liquids. The upper branch be- 
comes unstable to a FM state at kpa = 0.89(2), indepen- 
dent of the details of the interaction (in the zero-range 
limit). The critical kpa is similar for a purely repulsive 
interaction, but the result is non-universal and depends 
on details the potential; we will focus on hard spheres 
of diameter a. In both cases the critical value is consid- 
erably smaller than the Stoner MFT result (&fo)mft = 
7r/2 [3J. Despite similar values of the critical interac- 
tion, the behavior of the kinetic energy and the two- 
body correlations are qualitatively different for the upper 
branch and hard spheres. We also discuss the harmoni- 
cally trapped gas using the local density approximation 
(LDA). We conclude with a brief comparison of our re- 
sults with experiments [ij]. We find that some aspects of 
the experiment cannot be understood within our equilib- 
rium theory. 

Model: We consider a gas of N — (iV-f + A/jJ fermions 
of mass m with two species, denoted by "spin" f and I, 
which interact via a potential V(r). The Hamiltonian is 
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FIG. 1: (Color online) QMC energy per particle for the Fermi liquid state (for Nf = N± — 33 particles) as a function of fea 
for (a) the upper branch and for (b) hard spheres, compared with the perturbative result Eq. ©. Panel (b) shows QMC with 
and without backflow (BF) corrections. Backflow does not have a significant effect on the upper branch results in (a). There 
is a transition to a ferromagnetic state when the QMC energy crosses that of the fully polarized gas (horizontal lines), (c) The 
function r](r) describing BF correlation, discussed in the text, for the upper branch (UB) and hard sphere (HS) gas. 



with Tij — |r^ — rj-jj. For the QMC calculations we 
consider a cubic box with periodic boundary conditions. 
(At the end, we also discuss trap effects within LDA). We 
measure lengths in units of fc^ 1 , where kp = (37r 2 n) 1 / 3 
for a free Fermi gas of density n. We measure energies in 
units of €fg = 3£f/5 where ep = k F /2m (with h = 1). 

We consider two forms of the interaction potential. 
For the repulsive case, we use V(r) — Vq > for 
r < R and zero elsewhere. In the hard sphere limit 
Vq — > oo and the diameter of the sphere R = a 
scattering length. For the attractive case, we use 
V(r) — — (8 / mR 2 ) Vq /cosh 2 (2r I R) , extensively used in 
BCS-BEC crossover studies [ll|,ll2[. We choose the range 
such that kpR -C 0.1, so that we obtain universal results 
independent of the detailed form of V(r). We use Vq to 
tune the scattering length a > 0, such that V(r) has a 
single bound state. We will focus on scattering states to 
construct the upper branch wavefunction. 
QMC Results: The many-body wave function for a 
(paramagnetic) Fermi liquid is of the Jastrow-Slater form 
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The Slater determinants $fGct's are constructed from 
plane waves, and the symmetric Jastrow function f(r) 
accounts for interactions. 

We now argue that the upper branch Jastrow factor 
must be qualitatively different from the f(r) > used 
for the purely repulsive case. To see this, consider using 
a conventional nodeless /(r) for the attractive Fermi gas. 
This state is a normal Fermi liquid, and thus orthogo- 
nal (for large N) to the superfluid ground state 0, [ll[ 
of the BCS-BEC crossover for all kpa. However, the en- 
ergy per particle in this state is always lower than the 
free gas 36^/5, which means that the fermions do feel 
an attraction. In other words, this normal wave function 
necessarily has some pairing (bound state-like) correla- 
tions, and is therefore not on the upper branch. 



A necessary condition for a many-body state to be on 
the upper branch is that its energy per particle must 
be greater than 3ep/5. We must ensure that every pair 
of particles feels an effective repulsion. We achieve this 
by introducing a node in the Jastrow f(r). To deter- 
mine f(r), we use the lowest-order constrained varia- 
tional (LOCV) method (HI, which is well known in nu- 
clear physics and has also been used for strongly inter- 
acting quantum gases F2J, 1 14W . The LOCV equation has 
an upper-branch solution [15| f(r) with a node, whose 
location tracks the scattering length at small a [i.e., 
f(r) ~ (1 — a/r)] but then saturates at large a. 

We use QMC to calculate the energy for the up- 
per branch [Fig. HJa)] and the hard sphere Fermi gas 
[Fig. Hlb)] with Nf = N±. For small kpa, both results 
agree with the well-known perturbative result [2j 
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We note that Eq. ([3]) should be taken seriously only for 
kpa <S 1; the third order term is known to be non- 
universal, and depends on the detailed shape of the po- 
tential and on the p-wave scattering channel Q. 

A sufficient criterion for ferromagnetism (FM) is that 
the energy of the paramagnetic Fermi liquid state exceed 
that of the fully polarized state e^ G /(3e,p/5) = 2 2 / 3 ~ 
1.58. It is instructive to begin with simple analytical ap- 
proximations (even though these involve using Eq. ([3]) 
beyond its domain of validity!) The simplest approxima- 
tion is to just keep the first term in ([3]), the mean field 
Hartree shift. We find that this energy crosses that of 
the fully polarized e p FG at kpa = (2 2 / 3 - 1)9tt/10 ~ 1.66. 
This is slightly larger than the Stoner estimate of 7r/2, 
but still below the hard sphere solidification limit fcpa = 
(9-7T/4) 1 / 3 ~ 1.92. Including the second order term in (O 
increases the energy of the paramagnetic Fermi-liquid so- 
lution, and thus FM sets in closer to kpa ~ 1. 

The QMC energy for both the upper branch and hard 
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spheres implies a FM ground state for kpa > 0.9. We 
next address backflow to see how it affects our conclusion. 
Backflow: It is very important to include backflow 
which, as noted above, makes nontrivial modifications to 
the nodal surfaces and can lead to large quantitative ef- 
fects [§] in the ground state energy. Backflow modifies the 
single-particle plane wave orbitals </>k(i"io-) = exp [ik • Ti a ] 
used to construct the Slater determinants in Eq. ([2]) via 
the replacement r ia — > r icr + J2j 7 ?( r ij) r ij) where j labels 
particles of the opposite spin a. 

The optimal form of the backflow function ry(r) must 
be determined for each problem; it is known to be very 
different for the 4 He roton [l6| . for normal 3 He [l7j], and 
for the electron gas [ljj . Insight into the form of r\{r) for 
3 He came from analyzing the problem of a 3 He impurity 
in 4 He [llj- Following the same logic, we consider a single 
spin-down impurity in a spin-up Fermi sea. Omitting the 
details of our analysis (which will be reported elsewhere), 
we find the r/(r)'s shown in Fig.QJc). 

For the hard sphere gas the optimal rj(r) vanishes in- 
side the hard-core diameter and has a single peak jus t 
beyond it, qualitatively similar to the case of 3 He 0,0. 
As in 3 He, we approximate the form of r](r) by a Gaus- 
sian whose parameters we optimize. We use QMC to 
compute the energy of hard spheres using Eq. ^ with 
a nodeless Jastrow times backflow-corrected Slater de- 
terminants. We find that backflow leads to a significant 
lowering of energy [see Fig.[ljb)] that becomes more pro- 
nounced with increasing kp a. For example, there is a 
5.5% reduction in energy at kpa = 1. 

For the upper branch, we find that the form of the op- 
timal rj(r) is qualitatively different; see Fig. [lie). It is 
nonzero at the origin and decreases monotonically, with 
a power-law decay at large r. Further, T](r) changes 
very little with fcpa compared with the hard sphere case. 
The form of upper branch rj(r) is similar to systems 
with a soft-core, long-range repulsion, like the electron 
gas We use QMC to calculate the energy of the 

upper branch state (E|, with a Jastrow with a single 
node, times backflow-corrected Slater determinants. In 
this case the reduction in energy is small [Fig. QJa)] and 
falls within our statistical error of < 1%. 

We thus find that backflow is important for hard 
spheres when kp 1 is comparable to the hard-core diame- 
ter R = a. On the other hand, backflow effects are small 
for the upper branch, where kp 1 3> R, the range. 
Observables: For both the upper branch and for hard 
spheres, we conclude that ferromagnetism is energetically 
favorable, based on the crossing of energies of the para- 
magnetic Fermi liquid and the fully polarized FM; see 
Fig. [Ha, b). For the upper branch, we find that FM state 
is stable for kpa > 0.89(2). The order of the transi- 
tion requires a careful finite-size scaling analysis in the 
vicinity of the phase transition, beyond the scope of our 
present investigation. 

Although the total energies in the Fermi liquid phases 




FIG. 2: (Color online) Pair distribution function g(r) = 
g^l(r) for (a) the upper branch and (b) the hard sphere gas 
for various values of fcpa. Inset in panel (a) shows g(r — 0) 
as a function of kpa. 



in the upper branch and hard spheres are similar, the 
potential (V) and kinetic energy (K) are very different 
in the two cases. To understand this, it is illuminat- 
ing to look at the pair distribution function g-fi(r), de- 
noted by g(r) for simplicity. In the hard sphere case 
[Fig. Htb)], g{ r ) vanishes inside the hard-core and goes 
to unity at large separation. The potential energy (V) ~ 
J d 3 rg(r)V(r) then vanishes identically and the total en- 
ergy [Fig.QTb)] in the hard sphere case is entirely kinetic. 

In the upper branch, on the other hand, we find a large 
cancellation between a positive (K) and a negative (V). 
In marked contrast to hard spheres, the upper branch 
g{r) is extremely large at r — 0, has a pronounced dip 
at the node in the Jastrow f(r) and then goes to unity 
at large r; see Fig. [2][a). For the short-range attraction, 
the potential energy (V) ~ g(0) J d 3 rV(r) is thus large 
and negative, dominated by the growth of g(0) with in- 
creasing kpa [inset of Fig. Eta)]. This is compensated by 
a large positive kinetic energy (K) [Fig. EJa)] so that we 
find the total energy shown in Fig. [It a). 
Harmonic Trap: We first obtain from our QMC data 
the chemical potential (x = (dE/dN) as a function of 
density [Fig. Stb)]. We then invert this to find the equa- 
tion of state n{fi) of the homogeneous system. 

We restrict ourselves to the paramagnetic Fermi liquid 
regime here, and use the LDA [i(r) = //(O) — Vt rap (r) to 
study the effects of the harmonic trap Vt ra . p (r) with asso- 
ciated length scale clho- To compare with experiments, 
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FIG. 3: (Color online) (a) Kinetic energy (KE) per particle of 
the upper branch. Squares represent QMC data as a function 
k F a. The dashed line at fea = 0.89 shows the ferromagnetic 
transition at which the KE is greatly suppressed and then 
remains constant. The dot-dash curve is the LDA result for 
the KE versus k F a using the QMC equation of state. Within 
LDA, ferromagnetism appears at the center of the trap when 
k° F a ~ 1.1. (b) Chemical potential fx^kpa), related to the 
square of the LDA radius of the trapped cloud, for the upper 
branch. The perturbative result for [i = (dE/dN) obtained 
from Eq. © is also shown. 



we use the parameter k F — (24A r ) 1 ^ 6 /clho as a measure 
of the total number of particles N. To find the chemical 
potential at the center /i(0), we solve the LDA equa- 
tion (k F af = 2 3 / 2 48ir J^ 0) dft[fl(0) - /2] x / 2 n(/i). Here 
we have used dimensionless quantities jl = fi(0)ma 2 and 
n(p) = n(fj,)a 3 , where n(/i) is the QMC equation of state. 

We then find the density n(r — 0) at the center of 
the trap, from which we can determine the interaction 



parameter k F {0)a = [37r 2 n(/i(0))] 1/3 . We find that for 
k F a ~ 1.1, the trap center reaches kp(0)a = 0.89, the 
critical value in the homogeneous case. At this point the 
center of the trap should become unstable to ferromag- 
netism. We have also calculated within LDA the total 
and kinetic energies in the trapped system as functions 
of k F a; the latter is shown in Fig. |3ja). 
Comparison with experiments: While we were mo- 
tivated by the experiments of Ref . [4j , we focus only on 
"equilibrium" in the upper branch, and do not address 
dynamical questions. // three-body processes leading to 
molecule formation can be suppressed, there may be a 
window of time-scales where equilibrium physics in the 
upper branch, as described here, would be observed. The 
fcp a-dependence of g(r — 0) [inset of Fig.^a)] is relevant 
to the loss rate [ij} due to molecule formation. 

Even with these caveats, there are some aspects of the 
experiment which we can understand qualitatively and 
others we cannot. First, wc do find a ferromagnetic in- 
stability in the upper branch, but predict that it should 
happen in a homogeneous system at fcj?a = 0.89, which 
translates into the onset of FM at the center of the trap at 
k F a ~ 1.1, while the experiment sees interesting features 
only at k F a ~ 2. The behavior of the chemical potential 



saturates beyond the transition is qualitatively consis- 
tent with the experiment. However, the fcj?a-dependence 
of the kinetic energy is not; our results in Fig. OH a) are 
qualitatively different from the experiments. Finally, we 
have not addressed here the question of FM domains and 
their sizes, which is important to understand given that 
they have not been seen in the experiment. 
Conclusions: We show using QMC that fermions 
with effectively repulsive interactions become unstable 
to ferromagnetism beyond a critical interaction strength 
kpa ~ 0.9. This is true both for fermions in the upper 
branch (scattering state with positive a) of an attractive 
potential and also for hard sphere repulsion, despite im- 
portant differences in their short range correlations and 
the kinetic energy that reflect the underlying potentials. 
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Note added: As we were writing this paper, the work of 
Pilati et al. appeared [2(| • It addresses the same problem 
using a similar, but not identical, approach. Wherever 
they overlap, our results are in essential agreement. 
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